scalarField Rw (mesh.nCells(), 0.0);
scalarField Rn (mesh.nCells(), 0.0);
scalarField nRw (mesh.nCells(), 0.0);
scalarField nRn (mesh.nCells(), 0.0);
scalar CNVw=0, CNVn=0;
